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ABSTRACT 

The DEM L 316 system contains two shells, both with the characteristic sig- 
natnres of snpernova remnants (SNRs). We analyze Chandra and XMM-Newton 
data for DEM L 316, investigating its spatial and spectral X-ray featnres. Onr 
Chandra observations resolve the strnctnre of the northeastern SNR (Shell A) as 
a bright inner ring and a set of “arcs” snrronnded by fainter diffuse emission. The 
spectrum is well fit by a thermal plasma model with temperature ~ 1.4 keV; we 
do not hnd signihcant spectral differences for different regions of this SNR. The 
southwestern SNR (Shell B) exhibits an irregular X-ray outline, with a brighter 
interior ring of emission including a bright knot of emission. Overall the emission 
of the SNR is well described by a thermal plasma of temperature ~ 0.6 keV. The 
Bright Knot, however, is spectrally distinct from the rest of the SNR, requiring 
the addition of a high-energy spectral component consistent with a power-law 
spectrum of photon index 1.6-1.8. 

We confirm the findings of Nishiuchi et ah (2001) that the spectra of these 
shells are notably different, with Shell A requiring a high iron abundance for 
a good spectral fit, implying a Type la origin. We further explicitly compare 
abundance ratios to model predictions for Type la and Type II supernovae. The 
low ratios for Shell A (0/Fe of 1.5 and Ne/Fe of 0.2) and the high ratios for Shell 
B (0/Fe of 30-130 and Ne/Fe of 8-16) are consistent with Type la and Type 
II origins, respectively. The difference between the SNR progenitor types casts 
some doubt on the suggestion that these SNRs are interacting with one another. 

Subject headings: ISM: supernova remnants - ISM: individual (DEM L 316) - 
X-rays: ISM 
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1. Introduction 

The peanut-shaped nebula DEM L 316 was hrst noted by Mathewson & Clarke (1973) 
to have a high [S II]/Ha ratio, a typical signature of supernova remnants (SNRs). The 
authors designated the two lobes of this system shells A (the northeastern shell) and B (to 
the southwest); following this and subsequent work, we shall use the same designations. 
Mathewson et ah (1983) conhrmed that the DEM L 316 system had the optical, radio, and 
X-ray characteristics typical of SNRs. A multi-wavelength study by Williams et ah (1997) 
concluded that the most probable scenario for DEML316’s unique morphology was that it 
consisted of a pair of colliding SNRs. Evidence cited in support of this scenario included 
a change in radio polarization between the two shells, the kinematics of the two shells, 
multiwavelength morphological features and similar extinctions. 

Alternate scenarios for the physical association of the two shells include (a) that the 
shells are the result of a single explosion into a bi-lobed cavity, or (b) that the shells are 
isolated SNRs apparently juxtaposed along the line of sight, but in fact spatially separated. 
Nishiuchi et al. (2001), using a study of combined ROSAT and ASCA X-ray data, noted 
that the shells were spectrally distinct to an extent which would rule out the single-explosion 
scenario, and also concluded that, on the basis of its enhanced iron abundance. Shell A was 
probably the result of a Type la SN. 

In this paper, we present new Chandra observations of DEM L 316. These observations 
allow us, for the first time, to spatially resolve the structure of the DEM L 316 remnants, 
with reasonable spectral resolution over smaller regions. We supplement these data with 
XMM-Newton observations, which allow us to more reliably determine X-ray spectra for 
larger regions. §2 describes our observations, while §3 contains our analysis. A discussion of 
the implications of our Endings is given in §4, and those findings are summarized in §5. 


2. Observations 

We observed DEM L 316 with the C'handra Advanced CCD Imaging Spectrometer (ACIS) 
S3 back-illuminated chip (sequence number 500279, observation 2829, 40 ks). Data were 
reprocessed following procedures recommended by the Chandra X-ray Center (CXC): we 
removed the afterglow correction, generated a new bad pixel file, and applied corrections 
for charge-transfer inefficiency (CTI) and time-dependent gain for an instrument temper¬ 
ature of -120 C. We also applied background cleaning using the 5x5 pixel islands of our 
VFAINT mode observation. The dataset was filtered for high-background times and poor 
event grades, resulting in a total “good time” interval of 35.6 ks, and restricted to the energy 
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range of 0.3-8.0 keV, where the S3 chip is most sensitive. 

Spectral resnlts were extracted from the new event hie. Backgronnd regions were taken 
from source-free areas surrounding the SNRs, and the spectra of these background regions 
were scaled and subtracted from the source spectra. Individual spectra for regions of interest 
(Table 1), and the corresponding primary and auxiliary response hies, were extracted with 
the CIAO acisspec script and analyzed in XSPEC. Spectra were rebinned by spectral energy 
to achieve a signal-to-noise ratio of 3 in each bin. 

Regions within DEM L 316 were chosen to examine specihc morphological features seen 
in our X-ray images (Figure 1 and Table 1). Regions 1 and 2 enclose the approximate X-ray 
boundaries of Shell A and Shell B, respectively. Region 3 covers the bright central emission 
within Shell A, while Region 4 includes the limb of Shell A, excluding that central emission. 
Regions 5-7 cover brighter X-ray “arcs” to the north, east, and south of Shell A; Region 8 
samples emission between those arcs. In Shell B, Region 9 encloses the central emission from 
the remnant, and Region 10 covers a particularly bright knot of interior emission. Region 
11 incorporates the ring-like structure of brighter emission in Shell B, including the bright 
knot of Region 10. Region 12 covers emission on the limb of Shell B, excluding the center 
and the ring. 

To augment our spectral data for Shells A and B (Regions 1 and 2) we used XMM- 
Newton EPIC MOS and pn data from observation ID 0201030101, with exposure times 
of 11.2 ks (MOSl and MOS2 detectors) and 7.0 ks (pn detector). The pipeline-processed 
data obtained from the Science Operations Centre (SOC) were reduced using the Science 
Analysis Software package provided by the SOC. The data were hltered for poor event grades. 
Spectra for spatial regions corresponding to Regions 1 and 2 in the Chandra observations, and 
surrounding background regions, were extracted from the hltered event hies and rebinned 
to a minimum of 25 counts per bin. Joint hts were performed in xsPEC to the data from 
Chandra’s ACIS and the XMM-Newton instruments. 

For comparison, we utilized optical emission-line (Ha, [S II], [O III]) images taken with 
the Curt is-Schmidt telescope at Cerro Tololo Inter-American Observatory (CTIO), and radio 
images from the Australia Telescope Compact Array (ATCA) at 4.4 and 6.0 GHz (7 and 5 
cm). The angular resolution of the CTIO images was about 175; the half-power beamwidth 
of the radio images was 12". Other details of these observations are given in Williams et al. 
(1997). 
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3. Results 

3.1. X-ray Morphologies 

Previous X-ray analyses have noted a separation between the X-ray emission from the 
shells of DEM L 316. In the Chandra ACIS data, we indeed see two distinct shells; any 
overlap is too close to background levels to constitute a signihcant detection. Although 
diffuse emission between the shells is seen in the XMM-Newton observations, this emission is 
faint (~ 2a above the background with the EPIC pn); in addition, the SNRs are over 7' off- 
axis, and so the spatial resolution may be insufficient to distinguish betwen the boundaries 
of the two SNRs. 

Three-color images comparing the Chandra morphology to emission in Ha, [O III], 
and 7 cm radio (Williams et ah 1997) are shown in Figure 2a-c. In general, the brightest 
regions of X-ray emission appear to correspond to areas that are fainter at radio wavelengths. 
Conversely, the X-rays grow fainter toward the radio-bright limbs. Similarly, the X-ray 
emission appears well-conhned by the optical hlaments along the limbs of both shells. In 
Shell A, some hlaments across the face of the SNR appear to correspond roughly to regions 
of fainter X-ray emission, but there is no strong anticorrelation. 

Shell A appears to have a bright ring of inner emission, as well as bright “arcs” extending 
radially outward from that ring. As noted by Williams et al. (1997), the peak of the X- 
ray emission, now seen as the edge of this bright ring, corresponds to the location where 
optical hlaments (particularly well seen in [O III]) from Shell B overlap with Shell A. The 
southwestern side of the SNR is battened somewhat toward Shell B, and appears slightly 
brighter on the battened side. Notably, this battening occurs just north of the apparent 
overlap between the two shells, running parallel to the optical hlaments and radio-bright 
limb there. 

The faint, outermost emission from Shell B is irregular, and is elongated in the E-W 
direction. The primary source of this elongation seems to be an extension on the western 
side, which is matched by similar extensions in the optical and radio. The X-ray emission of 
Shell B is distributed over the face of the remnant, out to the boundary debned by the thick 
annulus of the radio structure. Corresponding to the inner radius of this radio annulus, there 
appears to be a brighter ring of X-ray emission, though this is not as prominent as that in 
Shell A. A bright X-ray knot is visible on the north side of this ring. Curiously, this knot is 
directly north of a “small-diameter source” noted in radio images by Williams et al. (1997). 

To examine the morphology as a function of energy, we created a three-color image (Fig¬ 
ure 2d) with emission divided into energy bands of 0.3-0.8 keV (soft), 0.8-1.5 keV (medium). 
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and 1.5-8.0 keV (hard). The middle band was selected to show the contribntion from the 
blend of Fe L lines that dominate the X-ray emission in this energy range, though it should 
be noted that other lines such as the helium-like blends of Ne and Mg also contribute to 
this band. The overall structure is similar in all three bands. Emission from Shell A shows 
considerably more medium-band emission than does Shell B; in particular, the central ring 
of Shell A is bright at these energies. In contrast, the X-ray Bright Knot in Shell B appears 
brightest in the hard band. A similar three-color map using the XMM-Newton data repro¬ 
duces these primary features at lower angular resolution: the medium band is much brighter 
in Shell A than in Shell B; and the location of the X-ray knot in Shell B is brighter in the 
hard band than is any other portion of either remnant. 


3.2. Spectral Fits 

We fit non-equilibrium ionization (NEI) plane-parallel shock models (“vpshock” in 
XSPEC)^, modified by photoelectric absorption, to determine the plasma parameters for 
each spatially-selected region. The model parameters are the mean shock temperature kT, 
ionization timescale (T=netion), elemental abundances (given as fractions of the solar abun¬ 
dance values of Anders & Grevesse 1989), and a normalization A proportional to the distance 
and volume emission measure (A(cm“^) = 10“^^ J nen]idV/A7rD‘^). Here, Ue is the electron 
density, uh the hydrogen density, tjon the shock age (time since the earliest material was 
shocked), V the volume occupied by the hot gas, and D the distance, all in cgs units. 


3.2.1. Overall characteristics of the two shells 

We extracted spectra from regions covering all of Shell A and all of Shell B (Table 1, 
Regions 1-2) from both the Chandra and XMM-Newton data. Spectra from the various 
instruments were fit jointly. The fitted spectral parameters are given in Table 2. It was not 
possible to obtain a statistically acceptable £t (at or above the 90% conhdence level) to both 
shell spectra with the same model parameters. 

The fitted absorption column densities were 3.6±0.6 x 10^^ H-atom cm“^ for Shell A, 
and 2.2 ±0.5 x 10^^ cm“^ for Shell B under a range of model combinations. Fits to Chandra 
data alone, however, gave absorption columns of 3.4±0.2 x 10^^ cm“^ for both remnants. 


^Details and references for the vpshock model can be found at 
http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/node40.html. 
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The value obtained from joint fits for Shell B is similar to that of 2.0-2.5 x 10^^ cm ^ 
estimated from hts to ROSAT data by Williams et ah (1997), although the value for Shell A 
is somewhat higher. The Shell A value, however, is reasonably consistent with the estimate 
of 3.8 xlO^^ cm ^ that those authors determined using H I data for the LMC (Rohlfs et al. 
1984) and a Galactic foreground of 5x10^'^ cm“^ (Schwering & Israel 1991). Data from a 
more recent H I survey of the LMC gives a higher hgure of 6.2x10^^ cm“^ (Staveley-Smith 
et al. 2003; Kim et al. 2003) toward DEM L 316, which when added to the aforementioned 
Galactic foreground gives a total absorption column of 6.7x10^^ cm“^. Note that the H I 
estimates include contributions from gas behind DEM L 316. Conversely, contributions to 
X-ray absorption from, e.g., molecular gas will not be accounted for in H I. 

Shell A is well £t by a single vpshock component with high iron abundance. Acceptable 
fits could be achieved with iron abundances between 1 and 5 times solar values; in all cases, 
substantially above the mean LMC abundance for iron. The Mg and Si He-like line blends 
are prominent in the spectrum, and fitted abundances for Mg, Si and S were well above 
typical values for the LMC ISM. The best-fit model yields an absorbed 0.3—8.0 keV flux of 
6±2x 10“^^ erg cm“^ s“^, an unabsorbed flux of 2±1 x 10“^^ erg cm“^ s“^, and a luminosity of 
3±2 X 10^^ erg s“^. The upper limit on an additional power-law component is an absorbed 
flux of 1x10“^^ erg cm“^ s“^, or 1.7% of the total flux; inclusion of such a component, 
however, does not improve the fit. 

For Shell B, a model with a single vpshock component provides only a poor fit. A 
substantial improvement in the fit is possible by using two components - either a combination 
of vpshock and power-law models or two vpshock models with the same abundances. For 
example, an F-test between the single vpshock model and a combined vpshock and power- 
law model gave only a 10“^^ probability of no improvement in the fit. The combination 
gives acceptable fits for values of the photon index T in the 1.6-1.8 range, within the typical 
range of photon indices for pulsar-wind nebulae in SNRs (Gotthelf 2003). Fits to Shell B 
do not require abundance fractions much above the LMC mean values. Models including an 
iron abundance of 1.0 solar or above could be ruled out at the 90% confidence level. The 
0.3-8 keV flux for this model is 8x10“^^ erg cm“^ s“^ (2x10“^^ erg cm“^ s“^ unabsorbed), 
implying a luminosity of 3x10^^ erg s“^. The power-law component takes up ~30% of the 
absorbed flux. 

The two-plasma fit gives an ionization-equilibrium component with a temperatnre ~0.6 
keV, and a non-equilibrium component with kT>3 keV. The htted abnndances are close 
to typical LMC abnndances. The fluxes and Inminosities for this model fit are the same, 
to one signihcant hgure, as those for the combined plasma/power-law fit. In this case the 
high-temperature component acconnts for ~50% of the hux. 



3.2.2. Regions within the SNR 


We used the Chandra ACIS data to obtain spectra from smaller spatial regions within 
DEML316 (Figure 1). We £t the best-fit spectral model found for Shell A (§3.2.1) to the 
various sub-regions within Shell A (Regions 3-8). For most regions (3-6,8) the emission did 
not significantly deviate from this model, and a joint fit (with only normalizations varied) 
produced statistically acceptable results. The exception was the S Arc region (Region 7), 
along the “flattened” side of the SNR. Emission from this region had an excess of soft 
emission when compared to the best-fit model for Shell A, and the joint fit to Shell A and 
the S Arc was poor. If we assume the abundances to be the same for this region as for Shell 
A generally, the S Arc region requires a slightly higher temperature (kT~1.5) and lower 
ionization parameter for a good fit. 

Similarly, we used the two-plasma fit to Shell B (§3.2.1) to perform joint fits with Shell B 
and its sub-regions (Regions 9-12). All of these regions were well described by the overall two- 
plasma-component Shell B fit, although the relative normalizations of the two components 
differ. The high temperature component accounted for > 70% of the unabsorbed flux for 
the Bright Knot (Region 10), ~40% of the flux for the Center Ring (Region 11) and < 30% 
of the flux in the Limb and Faint Center (Regions 9 and 12). 

We attempted to obtain additional information from these spatially resolved regions by 
fitting their spectra individually. To reduce the number of free parameters, we fixed the 
absorption column density to the mean value of 3.4x10^^ cm“^ found from fits to Chandra 
data for both shells, and fixed all plasma abundances except iron to 30% of their solar values. 
Given the low number of counts in most regions, it is unsurprising that the fits are highly 
uncertain, as seen from the broad errors on each fitted parameter. 

For all of these regions, a simple power-law fit alone could be ruled out at or above the 
90% confidence level. The best power-law fit, though still not statistically acceptable, was 
to the data from the Bright Knot region of Shell B (Table 1, Region 10). As with Shell B 
overall, an acceptable fit to the Bright Knot region could be obtained with a combination of 
plasma and power-law, or two plasma models at different temperatures. 

Broadly, the spectra of various regions within Shell A are quite similar. The tempera¬ 
tures are consistently within the range of 0.9-1.2 keV, and the fitted iron abundance within 
the range of 1.0-1.7 times solar values. The discrepancy between our fitted values for Fe 
abundance in the entire Shell A and sub-sections of this shell appears to be due to the dif¬ 
ference in fitting procedure. For Shell A itself, we allowed Mg, Si, S and Fe to vary, while 
for the sub-regions we only allowed Fe to vary. When we fix all abundances except Fe at 
0.3 for Shell A, we obtain Fe abundances of 1.7±0.2 solar. Allowing the S abundance, for 



example, to vary freely leads to a significant improvement of the fit (F-test probability for no 
improvement of 0.08) and a larger Fe abundance. The greatest differences between spectral 
hts to the various regions appear to be in the ionization timescale r, which is lower in regions 
along the remnant limb than toward the center. If we assume that most of the bright X-ray 
emitting gas is recently shocked, this would imply an increase in density toward the bright 
central region. 

For Shell B, the temperatures of the best-£t model combinations are generally in the 
0.6-0.8 keV range, with iron abundances less than 40% of solar values. The bright ring of 
emission (Region 11) is set apart from other areas of Shell B in requiring a high temperature 
or a second spectral component to provide an acceptable fit. The knot of emission embedded 
in this ring (Region 10) is also better described with the addition of a second spectral 
component, but the number of counts in this region is insufficient to distinguish between 
one- and two-component models at a high level of significance. 


4. Discussion 

4.1. Physical Properties of the SNRs 

We analyzed the overall emission from these SNRs, to determine abundance ratios be¬ 
tween key elements within the SNRs; compare our hndings with those previously made using 
ASCA and ROSAT instruments; and rehne our estimates of the physical properties of the 
SNRs based on the Chandra and XMM-Newton results. 


Abundances and SNR Types 

Our hndings conhrm the striking difference, hrst analyzed by Nishiuchi et al. (2001), 
between the spectra of the two shells of DEM L 316. Model hts to the spectrum for Shell 
A require an iron abundance signihcantly higher than that typical for the LMC, while no 
such condition is required for (or indeed permitted by) the spectrum of Shell B. The iron 
abundance found for Shell A, of three times solar abundances, is somewhat higher than the 
value of ~1.9 solar found by Nishiuchi et al. (2001). However, as discussed in §3.2.2, we can 
attribute this diherence to the dihering methods of spectral analysis. The hts of Nishiuchi 
et al. (2001) hxed all elemental abundances except Ne and/or Fe to the mean LMC values of 
0.3 solar. When we use the same procedure in htting the spectrum of Shell A, the htted Fe 
abundances of 1.7±0.2 solar are within the error range of the Nishiuchi et al. (2001) value. 
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The abundances found for Shell B are all well below the solar values, and in many cases 
near the typical LMC values (Russel & Dopita 1992). Our values for the iron abundance in 
the entire Shell B, and for sub-sections of this shell, all fall within the range of 0.16-0.4 solar. 
Nishiuchi et ah (2001) did not £t the Fe value directly for this shell, but noted that a fit with 
a fixed Fe abundance of 0.3 solar provided a reasonable spectral fit. These low abundance 
values suggest that the SNR is dominated by emission from swept-up material. However, 
the fitted values for O and Ne are sufficiently above their LMC values, even considering the 
large error ranges, to suggest that ejecta contributions to the spectrum remain significant 
for some elements. This case is far from unique, as ejecta enrichment is still detectable in 
other middle-aged LMC SNRs (e.g., Hendrick, Borkowski, & Reynolds 2003). 

The relative contribution of the lines from iron (compared with such elements as oxy¬ 
gen and neon) to the X-ray spectrum of a SNR is frequently used to discriminate between 
remnants resulting from Type la versus Type H supernovae (e.g., Hendrick, Borkowski, & 
Reynolds 2003; Nishiuchi et ah 2001; Hughes et ah 1995). Theoretical models (e.g., Badenes 
et al. 2003; Iwamoto et ah 1999) of Type la SNe and their remnants predict substantially 
higher ratios of Fe to other elements than those seen in remnants of Type H SNe (e.g., 
Iwamoto et al. 1999; Nomoto et al. 1997). A comparison of the ratios derived from Shells A 
and B to ratios from SN models is shown in Table 5. 

Using the abundances determined from the joint Chandra and XMM-Newton fits (Ta¬ 
ble 2) we find an overall 0/Fe ratio of ~1.5, and an Ne/Fe ratio of ~0.2 for Shell A. Models 
of Type la nucleosynthesis yields from Iwamoto et al. (1999) give typical 0/Fe ratios <1.0 
and Ne/Fe ratios <0.1. In contrast, modeled yields from Type H SNe from Nomoto et al. 
(1997), for progenitor masses of 15-40 Mq, give 0/Fe in the 8-370 range and Ne/Fe in 
the 0.6-30 range. Thus, the high iron abundance, particularly when contrasted with the 
less-enhanced Mg and Si, classes Shell A among LMC SNRs such as 0548-70.4 and 0534-69.9 
(Hendrick, Borkowski, & Reynolds 2003) and DEM L 71 (Hughes et al. 2003), which are ten¬ 
tatively categorized as remnants of Type la SNe due to high iron abundances. (Supporting 
this categorization, SNR 0548-70.4 and DEM L 71 also have Balmer-line dominated optical 
emission, a classic signature of a Type la remnant.) We therefore agree with the conclusion 
of Nishiuchi et al. (2001) that Shell A is probably a remnant of a Type la SN. 

For Shell B, depending on the model combination, the 0/Fe ratios are in the range 
30-130, and the Ne/Fe ratios in the range from 8-16. Typical ISM abundance ratios, based 
on the LMC values of Russel & Dopita (1992), would give 0/Fe of 13 and Ne/Fe of 2.4. The 
notably higher ratios found in Shell B suggest that these abundances not simply those one 
would expect from the swept-up interstellar gas, but are instead signihcantly enhanced by 
the contributions from ejecta. The abundance ratios for Shell B are clearly within the ranges 
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predicted for Type II SNe by the models of Nomoto et al. (1997). We therefore conclude 
that Shell B is very likely to be the remnant of a Type II SN. 


4 . 1 . 2 . Densities, Energies and Pressures 

The fitted parameters for the NEI plane-shock spectral models can be used to derive a 
number of physical properties for the SNRs. For this purpose we assume that hydrogen and 
helium are fully ionized, and that nue/nu = 0.1; thus, the total particle density is about 
1.92ne. We assume each SNR to be an ellipsoid with the line-of-sight dimension equal to the 
length of the minor axis, with volume V and hot gas volume hlling factor /hot- 

The spectrum of Shell A is reasonably consistent throughout the remnant, and can 
be well described by a single spectral component. We therefore use the overall spectral 
£t to Shell A as the basis for our estimation of its physical properties. For Shell B, the 
picture is complicated by the likely presence of more than one spectral component, and 
by variations over the face of the remnant. To characterize the properties of its overall 
expansion, therefore, we use the single-component spectral fit to the limb region (Region 12) 
as the basis for our calculations. From the normalization factor for a fitted spectrum, we 
can estimate the electron density and mass of the hot gas by: 

ne = 3.89 X cm-^ 


Mgas = 1.UnemnVfhot g 

The density and the plasma temperature (here given in keV) can be used to derive the 
thermal energy and pressure of the hot gas, according to: 

Eth = 4.60 X 10“®neTkevh^/hot erg 
Pth = 3.05 X 10“®neTkev dyn cm“^ 

Our derived values are summarized in Table 6, in terms of /hot where appropriate. 
Quoted errors are simply propagation of errors in fitted parameters and volumes, and do not 
take into account uncertainties due to the choice of spectral models and filling factor. 
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The properties dependent on Ug are affected by our choice of volume hlling factor. 
For spherical Sedov expansion, we expect /hot ~ 0.25 (Sedov 1959). Using this value, we 
hnd electron densities of 0.3±0.1 cm“^ for Shell A, and 0.28±0.03 cm“^for Shell B. These 
densities are consistent with the broad ranges found by Williams et ah (1997) from ROSAT 
data. The electron density found for Shell A is also consistent, within the errors, with the 
ion density (presumed equal to the electron density) found by Nishiuchi et ah (2001) from 
ASCA spectra, although it should be noted that their derived density was based on a value of 
/hot = 1 for this shell. If we similarly presumed /hot = 1 to reflect the hlled morphology, our 
density is only 0.2T0.1 cm“^ for Shell A. For Shell B, the value we hnd from the SNR limb 
is somewhat lower than that found by Nishiuchi et ah (2001); this difference is unsurprising, 
given that their spectral ht covered the entire remnant. 

The thermal energies for the hot gas, using /hot 0.25, are 1.5 ± 0.7 xl0^° erg and 
2.1 ± 0.2 xl0^° erg for Shells A and B, respectively. These values are also within the broad 
ranges found by Williams et al. (1997). The pressures found with this Filing factor are 1.4 
± 0.6 xl0“® dyn cm“^ for Shell A and 6.8 ± 0.7 xl0“^° dyn cm“^ for Shell B, somewhat 
higher than the ROSAT-derived values. We attribute this to the higher temperature (or high- 
temperature component) found from Chandra and XMM-Newton data; it is unsurprising that 
ROSAT, sensitive in the 0.2-2.4 keV range, would not have provided sufficient information 
at the high-energy end of the spectrum for accurate temperature determinations. (However, 
note that a higher filling factor for Shell A can lower the presumed pressure.) Our findings 
reenforce the conclusion of Williams et al. (1997) that the hot gas thermal pressure exceeds 
the sum of the magnetic pressure and the thermal pressure in the warm ionized gas seen in 
optical emission lines. 


4.2. Spatially Resolved Features of the Two SNRs 

4-2.1. Shell A 

Shell A fulfills many of the criteria for a “mixed morphology” (also known as “thermal 
composite”) SNR (Rho & Petre 1998): centrally brightened thermal X-ray emission, shell¬ 
like radio emission, and the absence of a prominent central compact source. The complex 
internal X-ray morphology does not preclude this interpretation, as structures similar to the 
inner ring of emission in Shell A have been seen in other mixed morphology SNRs, such as 
Kes 79 and HB 9 (Rho & Petre 1998). 

The identihcation of Shell A as a Type la SN might initially appear to contradict the 
mixed morphology interpretation, as these SNRs are usually associated with massive-star 



12 


phenomena. Mixed morphology SNRs are often found in the vicinity of molecular clouds, as 
seen in their association with maser activity (e.g., Yusef-Zadeh et al. 2003). Several mixed 
morphology SNRs also contain pulsars or PWN (e.g., Williams et al. 2005; Shelton, Kuntz, 
& Petre 2004; Rho et al. 1994). However, it is also not uncommon for SNRs from Type 
la SNe to be seen in proximity to massive-star phenomena. For example, the LMC SNR 
N103B, which shows X-ray inferred abundances suggestive of a type la origin^ (Lewis et 
al. 2003; Hughes et al. 1995), is found at the periphery of H II region DEM L 84, within 
40 pc of the NGC 1850 superbubble (Chu & Kennicutt 1988). Although not itself a mixed 
morphology SNR, N103B demonstrates that the association of many mixed morphology 
SNRs with massive-star phenomena does not rule out the possibility that such a SNR might 
have a Type la origin. 

Spectra of resolved regions for Shell A are remarkably similar in temperatures and 
abundances. This isothermal structure is similar to that seen in Galactic mixed morphology 
remnants such as Kes 79 (e.g.. Sun et al. 2004) and W44 (Shelton, Kuntz, & Petre 2004; Rho 
et al. 1994). The lack of large temperature gradients across SNRs are often attributed to 
thermal conduction and/or turbulent mixing (e.g.. Sun et al. 2004; Shelton, Kuntz, & Petre 
2004; Shelton et al. 1999). To see whether the former possibility is plausible, we examine 
the timescale for thermal conduction (Sarazin 1988; Spitzer 1962): 


^tc 


knj^ 


= 16,000 yr 


rZp 


/ 


1 cm“^ 7 \ 10 pc / \ 1 keV 


kT, 


- 5/2 


/luA 


where Kg is the Spitzer (1962) term for thermal conductivity (1.8 x 10“^ T^/^/lnA erg 
cm“^ s“^ K~^); and I the characteristic length scale for temperature variations (Te/|VTe|). 
The Goulomb logarithm In A is given by 


In A 


32.2 + In 



fcTe \ 

1 keVy 


which is ~33 for our overall Shell A values of Ue = 0.3 and Tp ~ T = 1.4 keV. We 
compare this estimate for the thermal conduction timescale (< 16,000 yr) to estimated 
ages for the SNR. Williams et al. (1997) suggest an age for Shell A of 27,000 yr, based 
on expansion velocities measured from echelle spectroscopy, and assuming Sedov expansion. 
Using the ionization timescale from X-ray spectral hts, Nishiuchi et al. (2001) estimated 
a larger age of ~39,000 yr. If the presence of magnetic helds does not signihcantly hinder 
thermal conduction, then, it is reasonable to expect temperature equilibration over the SNR. 


^This conclusion is disputed; see van der Heyden et al. (2002). 
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One of the primary factors that, it has been argued, may disrupt thermal conduction is 
the tangling of the SNR’s magnetic held by internal mixing. However, in that case, turbulent 
mixing (i.e. bulk motions in the gas) during the SNR’s development could itself produce a 
fairly hat temperature prohle (Cox et ah 1999; Shelton et ah 1999). In addition, one would 
expect these motions to create greater uniformity in the distribution of elements and in the 
ionization state of the gas. We would, therefore, still not expect dramatic spectral variations 
across the SNR if turbulent mixing is a prominent inhuence. 


4-2.2. Shell B 

In many respects. Shell B has the characteristics of an older remnant. The sub-solar 
metal abundances indicate substantial contributions from swept-up material to the X-ray 
emission, although, as discussed in §4.1.1, some ejecta contributions are still detectable. The 
low surface brightness of the X-ray emission and relatively slow expansion velocity are also 
consistent with this picture. 

As with Shell A, Shell B has some of the characteristics of mixed-morphology SNRs: 
it is generally shell-type in radio, lacks a bright outer limb in the X-ray regime, and has 
substantial internal emission. However, the case for Shell B is not as clear-cut as for Shell A. 
The surface brightness prohle of Shell B is shallower than that of Shell A, and the interior 
emission is more patchy. The presence of multiple spectral components for Shell B, and the 
indication that these components may be due at least in part to spectral variations across 
the remnant, are in marked contrast to the isothermal structure discussed above. 

The spatially-resolved spectra indicate that the high-energy spectral component is most 
prominent within the Bright Knot of X-ray emission. Given its proximity to the “small 
diameter radio source” mentioned in Williams et ah (1997), and the identihcation of Shell B 
as the result of a Type H SN, we must consider the possibility that the hard emission may 
come from a pulsar-wind nebula (PWN) embedded within the thermally-emitting shocked 
material. The high-temperature component of its X-ray spectrum is certainly consistent 
with the spectral properties of a PWN. However, the association of the X-ray knot with the 
radio source is uncertain, as the X-ray peak is offset from the radio peak by 13", greater 
than can be accounted for by the radio half-power beam width of 12". A search of the 
X-ray power spectrum of the XMM-Newton EPIC pn data for that region shows no evidence 
for periodicity in the emission; however, the low count rate of this source would make the 
detection of any pulsations very unlikely. 

In the absence of timing information or other indicators, and with the possibility re- 
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maining that the high-energy emission could be thermal, we cannot rule out the possibility 
that the Bright Knot emission may be simply a part of the remnant structure. For example, 
a dense clump of swept-up ISM could have been incorporated into the interior of the SNR. If 
this material were to evaporate into the hot cavity within the SNR, it could raise the surface 
brightness of high-temperature gas shocked during more energetic phases of the SNR’s ex¬ 
pansion (e.g.. White & Long 1991), and thus appear as a knot of higher-temperature X-ray 
emission. 


4.3. Implications for Collision Scenario 

Several scenarios have been put forward to explain the unusual bi-lobed morphology 
of DEM L 316. We can examine these scenarios in light of our new Endings. As put forth 
by Nishiuchi et ah (2001), we can decisively rule out the hypothesis that this system is the 
result of a single explosion into a bi-lobed cavity, on the basis of the wholly diEerent heavy 
element abundances found for the two shells. On the other hand, we certainly cannot rule 
out the possibility that the two SNRs are physically separate, and simply superposed along 
the line of sight. 

With regards to the central assertion of Williams et al. (1997) that DEM L 316 consists 
of a pair of colliding SNRs, the case for collision is weakened by the tentative conclusion 
that Shell A is the remnant of a Type la SN and Shell B is the result of a Type II SN. 
The explosions of two massive stars within the short lifetime of a SNR is conceivable if the 
two were members of a loose association, but if one SNR is the endpoint of the evolution of 
a low-mass star, no such association can be established. Since the observable lifetime of a 
SNR is relatively short (e.g., Slavin & Cox 1992), the probability of two independent SNe 
from non-associated progenitors within that lifetime is low. Even if the SN la progenitor 
had a main-sequence mass in the B-star range (Weidemann 1990), the required timescales 
for stellar evolution to culminate in a Type la versus a Type II SN are sufficiently different 
to make the probability of collision between their SNRs quite low. Thus, such a collision 
must be presumed to be a rare event, and the chances of observing it small. 

Interpretation of the X-ray morphology of the remnants remains ambiguous. Models of 
shock collisions often suggest that the hot cavities of two remnants undergoing a collision 
will be eventually connected by a “tunnel” allowing hot gas to flow from one to another, but 
this connection may not form until relatively late times compared to the SNR observable 
lifetimes (Bodenheimer et al. 1984; Jones et al. 1979; Ikeuchi 1978). Based on the clear 
spectral differences between the shells, we can eliminate the scenario whereby the hot cavities 
are interacting, although this leaves open the possibility of a collision observed before such a 
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tunnel has developed. In most cases the time between a SNR collision and the formation of 
a tunnel will be relatively short, so the liklihood of observing a system in this state would be 
small. However, if the SNRs collide during the late (radiative) stage of their evolution, such 
a tunnel may never form (Ikeuchi 1978). Such collisions would primarily be characterized 
by the dense “wall” of cool material along the interface between the SNRs. On the other 
hand, more recent laboratory and numerical simulations of colliding shocks (Velazquez et 
ah 2001) show that reflected shock waves from the collision of two shocks will push the 
cavities of hot shocked gas apart from one another, prodncing density structures reminiscent 
of the observed multiwavelength morphologies of DEM L 316. These reflected shocks would 
be expected initially to produce X-ray emission in the ring of dense gas along the interface; 
however, the very high density expected in this region also leads to rapid cooling of the gas, 
with a correspondingly rapid drop in the observable X-ray emission. In short, the Chandra 
observations do not provide any new snpport for the colliding-SNR scenario, but the scenario 
cannot be conclusively ruled out on the basis of these data. 


5. Summary 

Chandra ACIS images and spectra, snpplemented by XMM-Newton spectra, reveal new 
details of the DEM L 316 SNRs. 

1. Shell A exhibits a complex X-ray morphology, with a bright inner ring and “arcs” of 
brighter emission embedded in an elliptical region of diffuse emission. The brightest 
part of this ring occurs near an overlapping optically-emitting hlament from Shell B. 
The non-limb-brightened X-ray strncture is in marked contrast to the shell-like radio 
and optical morphologies, snggesting similarities to “mixed morphology” SNRs. As 
with many mixed morphology SNRs, little temperatnre variation is seen across Shell 
A. We suggest that the favored explanations for snch homogeneity in mixed morphology 
SNRs, thermal conduction and turbulent mixing, are also applicable to Shell A. 

2. Shell B also has its brightest emission in a ring-like structure well interior to the 
remnant limb. In particular, a small knot on this ring shows the brightest emission 
in the SNR. These regions are also spectrally distinct from the SNR limb, as they 
require the addition of a second, high-energy spectral component. The association of 
this high-energy component with the X-ray knot snggests the presence of an embedded 
pulsar-wind nebula, and its spectrnm is consistent with this snpposition, although we 
cannot rule out other possible explanations. 

3. We conhrm the Endings of Nishiuchi et al. (2001) that the spectrum of Shell A requires 
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a high iron abundance, and thus that this SNR likely resulted from a Type la SN. We 
are able to extend their work by explicitly fitting abundances for both shells, and 
comparing abundance ratios to those predicted by models of Type la and Type II SNe. 
We find that the 0/Fe and Ne/Fe ratios for Shell A are consistent with a Type la 
origin, while those for Shell B are consistent with a Type II origin. In the latter case, 
the observed ratios are significantly above those expected from swept-up ISM alone. 

4. The inferred physical properties of the hot gas are broadly typical of middle-aged 
SNRs. The relatively high thermal pressures in the hot gas for both SNRs emphasize 
the continued importance of the hot interiors to the evolution of the remnants. 

5. The observed spectral differences between the SNRs strengthen the argument of Nishi- 
uchi et ah (2001) that the shells are not part of a single bipolar SNR. We further 
observe that the spectrally-based SNR classifications weakens the case for the SNRs 
to be colliding, although the evidence is still far from conclusive. 


The authors thank the anonymous referee for comments which have helped to improve 
this paper. RMW acknowledges support from SAO grant GO2-3096A. This material is based 
on work supported by the National Aeronautics and Space Administration under NASA grant 
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Table 1. Spatial regions used for spectral fits 



SNR 

Region 

Counts'^ 

Center (J2000 RA, Dec) 

radii (") 

1 

Shell 

A 

entire SNR 

6352 ± 91 

05:47:21.4 -69:41:28 

63 X 49 

2 

Shell 

B 

entire SNR 

7085 ± 140 

05:46:58.8 -69:43:00 

104 X 64 

3 

Shell 

A 

Bright Center 

2822 ± 56 

05:47:19.0 -69:41:33 

21 

4 

Shell 

A 

Limb 

3024 ± 85 

05:47:21.4 -69:41:28 

63 X 49, -25'" 

5 

Shell 

A 

N Arc 

737 ± 31 

05:47:16.9 -69:40:58 

36 X 8 

6 

Shell 

A 

E Arc 

460 ± 24 

05:47:26.3 -69:41:16 

19 X 8 

7 

Shell 

A 

S Arc 

480 ± 26 

05:47:19.8 -69:42:10 

32 X 8 

8 

Shell 

A 

faint 

576 ± 31 

05:47:22.0 -69:40:58 

15 (N) 






05:47:24.2 -69:41:53 

15 (S) 

9 

Shell 

B 

Faint Center 

317 ± 22 

05:47:02.8 -69:42:55 

15 

10 

Shell 

B 

Bright Knot 

711 ± 29 

05:46:58.9 -69:42:27 

14 

11 

Shell 

B 

Center Ring 

3555 ± 79 

05:47:02.2 -69:42:54 

46, -16'" 

12 

Shell 

B 

Limb 

2714 ± 112 

05:46:58.8 -69:43:00 

104 X 64, -49^" 


^Background-subtracted counts in the 0.3-8 keV range, over 35585 s exposure time. 
'^Region with some interior emission excluded; radius of excluded region given as -X. 
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Table 2. Best-fit spectral models to SNR shells (joint XMM-Newton - Chandra fits) 


Parameter 

Shell A 


Shell B 


component 

vpshock 

vpshock 

Vpshock -h power-law 

vpshock -h vpshock 

Nff (cm“^) 

3.6lg;| X 10^1 

2.3j:o 2 X 10^1 

‘2.2t°Q% X lO^i 

2.lto j X 10^1 

kT (keV) 

1 4+0-3 

U.DO_q Q2 

u.oy_Q Q 3 

0.57l°;j 51® 

0/00 

0.22t°;® 

Q y+0.2 
0'-'-0.2 

0 75+0 07 
'-'•'O_0.07 

0.4l°:j 

Ne/Ne 0 


0.6tg;i 

o+DI 

0.7tg;i 

Mg/Mg 0 



o.elgi 

n 7+0-3 

Si/Si 0 

0 Q+0-4 

0.351H 

U.O_0.2 

'^•0-0.1 

Fe/Fe 0 

r\ -T-I-O.S 

- 0.5 

0 10+^'33 
h'--LU_Q Q2 

n 1 c+0.05 

0 .iO_Q Q 3 

rv (^o-t-0.09 

r (cm”^ s) 

l-+a5 ^ 10“ 

< 5 X 1013 

<5x1013 

< 5 X 1013 l-StoJ ^ loll 

r 



17+0-4 

-*-+-0.2 


norm (cm“^) 

X?ed 

dof 

2.lto;? X 10“"' 
1.19 

210 

l-+al X 10-3 
1.36 

507 

1-411o o 4 X 10-3 5.3lo;| X 10-5 

1.20 

505 

S.lto s X 10-1 S.ltjJJ X 10-1 

1.17 

504 


Note. — Spectra cover the range between 0.3-8.0 keV. Elements not listed are fixed to a LMC mean abundance of 30% solar 
(Russel Sz Dopita 1992). Quoted errors are the statistical errors in each fit parameter at the 90% uncertainty level. 


Table 3. Best-fit spectral models to regions in Shell A 


Parameter 

Brt Ctr 

Limb 

N Arc 

S Arc 

E Arc 

Faint 

kT (keV) 

0 Q1+0-01 

-1 90-1-0.07 

Q 5 

0.9812:22 

4 - 4 _ 0.1 

1 + 2 
-^-0.2 

0 Q 1 +0-09 

o-yi_o.oi 

Fe/Fe 0 


1 4+0-4 
-*-•^-0.2 

-| 9-1-O.5 
-*-•^-0.4 

i+O.S 

-*^•^-0.3 

1+2 

-*^-0.5 

1 4-1-0.9 

4-^-0.4 

r (cm”^ s) 

4lL X 10^" 

2.2j:2'2 X iQii 

112°4 X I0I2 

2C1 X loll 

6lf 1 X lOii 

5I2 X loll 

norm (cm~^) 

X?ed 

dof 

2.091217 X 10-1 
1.19 

91 

2.ol2;j X 10-1 

1.17 

177 

6.412:22 X 10-5 

1.04 

54 

3.812:21 X10-5 
1.34 

44 

3.612:2 X 10-5 

1.03 

36 

4.312:2 X10- 

1.27 

57 


Note. — Here and in Table 4, Nff is fixed at 3.4x10^^ cm ^ and non-Fe abundances at 0.3 solar. The “Limb” region 
incorporates the “Arc” and “Faint” regions. All model fits in this table use the vpshock model. 














22 


Table 4. Best-fit spectral models to regions in Shell B 


Parameter Limb 


Bright Ring 


Faint Ctr 


Bright Knot 


component 

vpshock 

vpshock 

vpshock + pwrlw 

vpshock 

Vpshock 

vpshock + pwrlw 

kT (keV) 

0-78tH5 

-*-•^-0.2 

0 53+0 04 
h'-oo_o 03 

0 

0 43+0.09 

n 65+0-08 

‘J-D^>_0.09 

Fe/Fe 0 

0.29in? 

0 Qi +0 04 

h'-Oi_o.o5 

0 16+003 

h'-J-D_o 04 

0 

^•4-0.1 

'^•^-0.1 

0.3^ 

r (cm”^ s) 


sli X ioi° 

6lf X 10“ 

< 4 X 10i5 

< 4 X 1045 

< 5 X 1045 

r 



9 9+0.4 

z.z_o 2 



^•’^-0.4 

norm (cm~^) 

d.eto J 10“^ 

X 10-^ 

8.3tJ^l X 10-^ 

l.ll°l X 10-4 

5.+° j X 10-5 

etf X10-5 

norm2 (cm~^) 



6^1 X10-5 



2.91^1 X 10"® 

X?ed 

1.22 

1.26 

1.05 

0.80 

1.29 

1.08 

dof 

268 

189 

187 

35 

61 

60 


^Could not be adequately constrained, so fixed to LMC mean abundance. 


Table 5. Element ratios for SNR shells 


Ratio 

Shell A 

Shell B 

LMC'* 

Type la models'^ 


Type II models^ 






W7 

WDDl 

15M0 

2OM0 

25M0 

4OM0 

O/Fe 

1.5 

30-130 

13.2 

0.67 

0.46 

8.48 

66.0 

178 

367 

Ne/Fe 

0.2 

8-16 

2.4 

0.016 

0.0060 

0.59 

9.08 

29.9 

23.0 

O/Si 

5.9 

19-48 

3.5 

1.63 

0.57 

8.18 

25.8 

45.1 

30.6 


^Typical LMC abundance ratios (Russel Sz Dopita 1992) 

*^Based on nucleosynthesis models for Type la SNe (Iwamoto et al. 1999). W7=fast defiagration model, WDDl=slow 
defiagration, delayed detonation model. 

^Based on nucleosynthesis models for Type II SNe (Nomoto et al. 1997). 


Table 6. SNR properties from spectral hts 



Shell A 

Shell B 

He (cm ^) 

0.16 ± 0.07 4 -y" 

0.14 ± 0.01 4 -y" 

M^gas (^©) 

50 ± 20 / 4 /f 

no ± 10 

Eth (erg) 

3.0 ± 0.9 xl050 

4.2 ± 0.5 xl050 

■Pth (dyn cm-2) 

6.8 ± 0.9 xl0-i° 

3.4 ± 0.4 xl0-i° 
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Fig. 1.— Adaptavely smoothed Chandra ACIS image with spectral regions marked and 
labeled as listed in Table 1. Regions with dotted lines are excluded from larger regions. 
Excluded point sources are also shown with dotted lines. 
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Fig. 2.— (a) 3-band image showing Ha from the CTIO Curtis-Schmidt telescope (red), 7 cm 
radio map from the ATCA (green), and adaptively smoothed Chandra ACTS image (bine), 
(b) 3-band image showing Ha (red), [O III] (green), and smoothed ACTS image (blue), (c) 
2-band image showing ATCA 7 cm (red) and smoothed ACTS image (blue), (d) 3-color 
Chandra ACIS image showing emission at 0.3-0.8 keV (red), 0.8-1.5 keV (green), and 1.5-8.0 
keV (blue). The X-ray images were adaptively smoothed on the same scale. 











